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ABSTRACT. Maize Lethal Necrosis disease (MLND) is a viral disease that can cause fatal 
damage to the crop of maize plants. This is very common in East Africa countries and 
Democratic Republic of Congo (DRC). In this manuscript, a mathematical model has been 
developed to study and analyze the dynamics of the MLND in the maize crop population. 
The disease free (DFE) and endemic equilibrium (EE) points of the model has been com- 
puted and the basic reproduction number (Ro) derived using the next generation matrix 
method. We performed sensitivity analysis by using parametric values from literature and 
estimated ones. We found that the rates of transmission, Amo, Go and mm are the most 
positively sensitive parameters. Numerical simulations were also performed to verify the an- 
alytical results.Thus, this research work recommends that deliberate strategic intervention 
should be targeted on the disease transmission rates which are significant for MLND trans- 
mission in order to eradicate the disease or reduce the intensity of the disease transmission 
in the maize population. 


1. INTRODUCTION 


Research teams from various countries in Africa are under pressure of finding solution to 
threaten food security which is due maize lethal necrosis diseases (MLND). MLND is a new 
serious disease which emerged for the first time in Kenya in 2011 and spread to Tanza- 
nia and Uganda in the year 2012 [1]. In these occasions a number of damage and losses 
were recorded. The disease is caused by maize chlorotic mottle virus (MCMV) and sugarcane 
mosaic virus (SCMV) [2] or co infections of MCMV and any other types of viruses such 
as Wheat Streak Mosaic Virus (WSMV) or Maize Dwarf Mosaic Virus (MDMV) [2] form 
MLND. 

The transmission of MLND is attributable to various insects such as maize thrips, rootworms, 
leaf beetles and leaf hoppers. However, research findings indicate the disease transmission 
through seeds normally occurs at very low rates [2]. Also, the MLND virus may spread 
through soil, infected plant debris and aphid vectors [3, 4, 5]. 

Mathematical modeling has also been a significant mainstay in the development of strategic 
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intervention programs towards disease control in the plant crops population. Several deter- 
ministic models have been established to describe and analyze the dynamics of diseases in 
plant populations,but few correlate to this studies. Thus mathematical modeling has not fully 
been employed as an alternative way to describe and analyze the dynamics of the disease 
under consideration in our case. We present a deterministic model to study and analyze the 
dynamics of MNLD in the maize population in this paper. We believe that the results of our 
research work will be useful in finding suitable means of controlling the disease transmission, 
or rather eradicate it.This may ensure farmers and peasants maximum maize harvest for 


food security. 


2. MODEL FORMULATION 


This paper presents an SI-SEI-type model of host and vector populations that incorpo- 
rates: The maize host population, which is categorized as Susceptible sub-population (Sm); 
Exposed maize sub-population (Em) and Infected maize sub-population (Im). The model in- 
cludes: the vector population, which is categorized as Susceptible vector sub-population (S,) 
and Infected vector sub-population (J,). It also includes: virus in the environment (P,). In 
the model, there is no recruitment rate since it describes a single season and no plantation of 
maize plants. The class Sm declines constantly at the rate (7), which is due to force of infec- 
tion between virus in the environment (P,); direct contact of maize to maize (Im) and vector 
from infected maize plants (I„). The infected maize plants finally increases the number of 
exposed plants (Em) at the same constant rate (7). However, some susceptible maize plants 
remain in the exposed state while the remaining ones immediately progress to infected maize 
class, (Im) at a constant rate (am). Then the infected maize has disease-induced constant 
death rate o,,.There is no natural death in the host plants population for all three classes as 
maize plants never naturally die till harvested at the end of season. The infected maize and 
the environment transmit the MNLD virus to the susceptible vector, S, at constant rates 
(Gmv) and (8o) respectively.The susceptible vectors are constantly recruited at the rate A; 
constantly progress to infected vector class with mass action A(P,, Im), with proportionality 
constant rate À and die naturally at a constant rate u. The infected vectors also naturally 
die at the same constant ratej:,. The MNLD virus in the environment are recruited at the 
rate proportional to (@P,) and from the infected maize plants at a constant rate (o), and 
then die naturally at the rate( uo). 

In the formulation of our model the following assumptions were considered: 


i Pathogen from infected maize plants reaches the vegetative environment by vectors 
and through shedding which is due to wind, people, rain and birds. 


ii Transmission of infection from the vegetative environment to susceptible maize is 


through vectors and shedding which is due to wind, people, rain and birds. 
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iii Some maize plants are exposed but resist the MLND throughout the season. 
iv Once the vectors become the carriers of pathogens, it is for the whole of their life. 


v The population is heterogeneous. That is to say, the individuals which compose the 
population can be grouped into different classes in accordance to their epidemiological 
state. 


vi Each susceptible individual in the class has equal chance to be infected by contagious 
individuals if it happens they come in contact. 


vii There are no immigrants and expatriates (emigrants). The only pathway of entering 
into the population is through sowing seeds and the only way of exit is through death 
from MLND-related causes. 


2.1. Model variables and parameters. A complete explanation of the variables and pa- 
rameters that are used in the model are summarized and described in Tables 1 and 2, 




















respectively. 
Table 1: Variables and their description 
Variable Description Units 
Nm Total number of maize plants population plant ha~ 
Ny Total number of vector population vector ha! 
oe Susceptible population maize plants plant ha~ 
Sy Susceptible vector population vector ha~! 
T Infected population maize plants plant ha~* 
Iy Infected vector population for maize plants vector ha~! 
Em Exposed maize plants with latent incubation period plant ha~ 
P, Other plants carrier of pathogens (vegetative environment) cells ha~! 
Table 2: Parameters and their description 
Parameter Description Units 
Bris Contact rate from infected maize to susceptible vector 
population ha cell! day} 
lo Natural death rate of vector within environment ha plants! day“! 
Be Contact rate of other plants (not maize) carrying 
pathogen with maize ha cell! day} 
Bov Contact rate of virus in the environment to the 





Continued on next page 
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Table 2 — Continued from previous page 








Parameter Description Units 

susceptible vectors ha cell~! day~* 
Diri Contact rate from virus carrying vectors to Sm maize 

population ha cell-! day™t 
Dmm Contact rate of contagious maize to susceptible maize 

population ha cell~! day™t 
Avo Contact rate from virus carrying vectors to the vegetative 

environment ha cell~! day“! 
Amo Contact rate from infectious maize plants to the vegetative 

environment ha cell~! day~* 
A Recruitment rate into susceptible vector 

population cells day~* 
n Force of infection rate of 3 masses of action due 

vectors, Im and virus in P, ha cell-1 day™t 
À Force of infection rate of combined masses of action 

due infected maize and virus in the environment. ha cell~! day~* 
Le Natural death rate of vector population ha plants! day! 
Om Death rate due maize infection plants™! day™t 
o Contribution of infected maize to the growth of 

viruses. cells plants! day~! 
0 Contribution of infected plants to the growth of viruses 

in the environment cells plants! day~! 
Am Progression rate of exposed maize plant to be 


infectious plant day™!t 
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2.2. Model flow chart. Using assumptions of the model and defined 
variables/parameters, the dynamics of maize lethal necrosis disease (MLND) in maize 
population can be shown as the flow chart in Figure 1 depicts. 
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FIGURE 1. Schematic of MLND with susceptible maize plants population 
which interacts with forces of infection from viruses in the environment (P,), 
infected maize (Im) and infected vector (/,,) populations. Epidemic starts from 
susceptible maize (Sm) and moves to exposed class (Em) and infected class (Zm) 
respectively. 


2.3. Equations of the model. Considering the compartmental diagram as described by 
Figure 1, we formulate basic mathematical model which shows transmission dynamics for 
MLND using the following differential equations: 





(1) TE = (Ber Ep, + Bom se + Bam E) Sm: 

2) TP = (Bagh + Bom wR |S — Om Em 
(3) “mn SOF y= Gila 

(4) We = A (Baz er 25 + Bm RE) Sv — Ho 

(6) TE = (Bag ep + Pm E) S = polo 

(6) TE = Amott + Ag = HoPo 


with initial conditions: 
Sm > 0, Em 2 0, Im 2 0, P, > 0, So > 0, Im 2 0. 
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Note that the total number of maize plants population is given by 
(7) Nm = Sm + Em + Im, 

while the total number of vector population is given by 


(8) N, =S, +L. 


2.4. Basic properties of the model. In order to see if the model is mathematically and 
epidemiologically well posed to study the MLND dynamics, we considered two basic proper- 
ties of the model. We used Box invariant of Metzler Matrix to show the existence of invariant 
region |6, 7] and basic standard method to prove the positivity of solutions. 


2.4.1. Invariant region. We proof the invariant region as follows: The system of the model 
(1) can be written as 

s = A(X) X +F with X=(Sm, Em, Im, Sv, Iv, Po)? and the constant term F = (0, 0, 0, ,0,0)*. 
That is 








=F 0 0 0 0 0 

Fm Am 0 0 0 0 
0 m —Om 0 0 0 

A(Xy= ý i 
0 0 0 —Fy — fly 0 0 
0 0 0 Fy Ly 0 
0 Q Ame 0 ie 
and 
ie, ee a 
Fr SRT P i Pom ay i Pnom) 





Pe aa oles 
Wa = (Bove. 1 Brox): 


A(X) is a Metzler Matrix Y Xe RÊ which all off diagonal terms are non negative and 

F > 0.The system “4"=A(X)X +F is positive invariant in R. We can conclude that the 
feasible region Q is a set of Q=(Sm, Em, Im, Sv, Iv, Po) € Re with initial data 

Sm > 0,Em > 0,Im > 0,5, > 0,» > 0,P, > 0. Hence, the solution stays in the region if it 
started in the region. 

Similarly, letting the initial data be {Sm (0), Em(0), Im(0), 7,(0), P,(0)} € Q and using stan- 
dard method, the solution set { Sm, Em, Im, Sv, Iv, Po} of the model system (1) was shown to 

be non negative V t > 0. It was therefore concluded that the solution set Q={ Sm, Em, Im, Sv, Iv, Po} 
of the model system (1) is non-negative for all t > 0 and so it is sufficient to study the dy- 
namics of MLND. 
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3. MODEL ANALYSIS 


3.1. Existence of disease free equilibrium point. The disease free equilibrium (DFE) 
point of the model system (1) is given by 


A 
(9) EP = (Sm, Em, Im, Sv, Io, Po) = (Nm, 0,0, —, 0, 0). 
When there is no disease, there will be susceptible maize population plants, Nm. On the 
other hand in the susceptible vector population there will be new reproduction of vectors A 
and dying Hv, 
where Nm is the constant with adoption value of 44,000 [8] maize population in one hector. 


3.2. The basic reproduction number R,. Analysis of the equilibrium point, is carried 
out by considering spectral radius of the matrix FV~! called threshold quantity. The basic 
reproduction number is the measure of secondary contagious number as a results of one 
contagious maize plant in the susceptible maize population. 

With the help of Next Generation Matrix, the approach adopted, we determine the basic 
reproduction number [9]. We consider the population distinguished by six different classes 
which are Sm, Em, Im, Sv, Iy and P}. 

Firstly we arrange the system to get group of infectious classes only that is (Em, Im, Iv, Po). 
Once again we assume f;(x) be the rate of introduction of new infectious (transmission) in 
compartment i,v; (2) be the transmission after new infectious (transition rate by all other 
means) and v; (x) rate of transfer of individual out of compartment j. 

The disease transmission model comprised of the system of equations 


X; = f(z) — v;(2), 
where 
vj(@) = uf (x) — v} (2). 
To obtain the matrices F and V of dimension ’n’by’n’ we differentiate vectors fj(x) and 
vu;(x) respectively. 








_ OF; (x) 

_ v(x) 
with 

I<, i<n, 


where xo is a vector column from DFE. From the model equations 
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En Ob Tin 
(filo) = |p | ands e me? , 
Im Iv 
0 Amon + Aov Ne — oP, 





Ee Ge f biar H Bren Sm 
Fy = Neeson Brune 


) 
Differentiate (Fm, F,,0,0)" with respect to (Em,Im,I,,P,)? and obtain 





m m k m 
0 Cie Bian Bone Fn 


0 0 0 0 
F= 0 By, Se A P 
MU Nym ov (k+P)2 
0 0 0 0 
At disease free Sm = Nm which is the number of maize grains sowed and germinated in a 


maize plant, and S, = a as from 


ay =A-— (Baar, } Pm dN = ore 


Hence matrix F reduces to 








0 pmm m Ben 
0 0 0 0 
F= 
0 Bain 0 Baier 
0 0 0 0 
Am Em 
Again differentiate v;(x,) = k - a 
Ly v 


Niles al Nein Fe Pals 
with respect to (Em, Im, Iv, P)” to obtain matrix V. 
Hence, 





tees —Am Om 0 0 

0 0 [by 0 

OS Fae ae A 

We find the inverse of V i.e 
E 0 0 0 
eae: + m 0 0 
7 0 0 ex, . i) 
Hv 

Amo Amo Avo 1 





2 
3 
Q 
3 
= 
ic} 
2 
3 
Q 
3 
T 
o 
Z 
B= 
e 
cs 
or 
= 
= 


MLND DYNAMICS AND SENSITIVITY ANALYSIS IN A MAIZE POPULATION 9 


And multiplying matrix F and V~! we obtain 











bmm -+ BoAmo Bmm | BoAmo BumNm | AvoBoNm BoNm 
Om komo Om ' kombo Nvo N kNoHvho kuo 
ae 0 0 0 0 
Bmv A | AmoBovA BmuA | AmoBouA AvoBovA Bov A 
NmOmby kNmOmbv Lo NmoOmbv kNmOmbv Lo kK Ny u2 Mo khv Ho 
0 0 0 0 


According to [10] the matrix can be reduced into 2 by 2 matrix. [11] The dominant eigen- 
value of this matrix is Ro, which can be obtained from the trace and determinant of that 














matrix as 
Bmm J Bomo Bum Nm | BoNm>2vo 
K = Om ! kom Ho Noho l kKNoy Hv Ho 
Bmw A | Bov A Amo Bov A Avo 
NmHvam ! klv Nmam Ho ko Nv Hv? 
pmm A Gone AvobovA 
Trace(K) = | | 5 
Om kompo =k Nope to 


or in a simple form: 


NA Ce alts + Bomo) + Avobov A 








T K= 
race(K) EN iP ii, 
A kbmubumho = Dri ouawe + BmuboAvo + Pov bomAmò 
Det(K) = — ( Noe ) 


But according to [10] we have this formula for 





R, = p(K) = ; (Trace(K) + /Trace(K)? — 4det(K)) 


We substitute trace and determinant in the formula above and obtain a condensed form as 


shown below. 











1 mm ee A Nvo ov 
R, = p(K) = p | p | 
2 \ om  KOmto N, kĉžho 
i 1 Bmm BoAmo i A AvoBov 2 i A kBmv Bum Ho a BmmPBovrAvo + BmvBoAvo + BovBumAmo 
(to) | s( om ij komo Ny Ere) l 4N, ( kohv? am ) 


3.3. Endemic equilibrium point. Endemic equilibrium point of the model system (1) 
E* is a steady state solution where by the disease persists in the population or an equi- 
librium comprised of diseased sub-populations i.e. Em 4 0, Im # 0, Ty # 0 and P, F 0. 
By equating all expressions for the equations in the model system (1) to zero and let E* 
=(S* , Ek, IZ, S*,1*, P*) we compute the system and arrive to the cubic polynomial equa- 


tion as follows. 


(11) P(I*) = Ar + BI? — DI* -C 


where 
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A= Bram loNo mel Ga + om), 

B = (Am + Om) boh Ho NE NmAmo + BmvbtultoNmAmoA(Om + Om) + BmmbatoNe Nmk(Om + Om) 
+ Brmbtobo NmNyXovA(Am + Om) > —BmmbetoNmNyAmom — Cinti N, AmXous 

D = —Bobolty NZ, No AnA (Gm + Om) + Bmoth NaN AOmk (Gm + Om) 

= Prutl NON Onk An Fom) = Bo HoN No Amoi 

— Bmmlivtio N? N, Aam — bmnp UNZ N? kam 

— poptoN p NyAmoAGmOm — hoho NaN komam, 

C= e Boao N n NyA N a — Briha N n No Akm — Orita N n NyA A On 


The obtained polynomial function with cubic degree indicates the existence of endemic equi- 
librium points. Using Maple 18 (32 bit) we solve the polynomial equation to obtain three 
distinct roots of which one is real and the other two are complex. 

We rely on the real root and define an equilibrium point which is comprised of epidemic 
populations and given by the endemic equilibrium point as 























S5 Ay 
E* Ap 
e I, Ag 
(12) (E°) = f= 
oe Ag 
I As 
P Ag 
where 
Am + om af? — bag — 6a? fl 
iaia ee 
Am 6a? f 
Om (af? —6ag — 6a? fl 
No i ; 
Qm 6a? f 
Af? — 6ag — 6a? fl 
As = , 
6a? f 
A 
A, = 
Bov P: I Bmvlin | : 
k+P, | New Hv 
2 
Bov Po | Bmolh | 
(gek PAN Ht) 
A5 m ; 
Apy 
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and 











1 4/108 CA? +36 (D) BA + 12 V3\/27 A2C? + 18 ABC (D) + 4 AD? — 4 B?C — B? D? A — 8 B? 














6 A f 
g= A108 C42 + 36 (D) BA + 12V3,/27 A?C? + 18 ABC (D) 4 AD’ — 4 B3C — B2 D? A — 8 B3, 
2 
a= A. 


3.4. Sensitivity analysis. Sensitivity analysis can applied to test the ability of the model 
in relation to prediction of model parameters[11]. In this section we use sensitivity analysis 
to determine the impact of parameters on R,. In order to determine an excellent way that 
can reduce maize death and morbidity due MLND, it is crucially important to understand 
the proportional importance of factors that are liable for the transmission and prevalence of 
the infection. 

Table 3 shows the parameter values of maize disease model. The parameters are taken 
from the previous studies that correspond to this study, existing information and through 
approximation. 

To perform sensitivity analysis of R, we need to first determine sensitivity indices’s for each 
of parameters. 

We define sensitivity indices’s as partial derivative of R, with respect to a given parameter 
as multiplied by the ratio of given parameter to the basic reproduction number Ro. 
Mathematically, this definition is given by: 

OR, Mi 

—)x —. 

Om; Ro 


The approach used in equation (13) above is called normalized forward sensitivity indices of 


(13) rie 


Mi 


R, with respect to parameters m;, where m; is any parameter from Ro. 


3.5. Parameters estimation and adoption. In this subsection, we will review some jour- 
nals to get proper parameters and where necessary approach of estimation will be used. 

The primary infection rate of other plants population rather than maize can be greater or 
equivalent to 0.02 i.e 0.04,0.06,0.08... [12] 

Transmission of Maize Chlorotic Mottle Virus by Chrysomelid Beetles and its cleoptera 
species of 8 classes was estimated [13]. From acquisition to inoculation is 8 days. The rates 
of these different classes were 0.390.0.905,0.052,0.024,0.170,0.139,0.156 and 0.Finding the av- 
erage we obtain 0.2295 and again divide by 8 days and results to our one of transmission 
rate Bmm 0.0286875. 

Leafhopper transmission rate of ACMV 0.020(standard) 0.005 to 0.100 (range) (diseasedplant)~‘day7' 
and transmission rate of EACMV 0.030 (standard) 0.01 to 0140 (range) (diseasedplant)~‘day7' 
[14]. Although MCMV could be transmitted mechanically and by seeds at very low rate 


12 ALOYCE, KUZNETSOV AND LUBOOBI 


(0.008 0.04), the insects have crucial role in disease transmission and as the matter of fact 
leafhoppers are responsible by settling on the crops{15]. Therefore, from this description we 
can estimate some of our model parameters between such interval as from 0.005 to 0.100 as 
we can see in the Table 3. 

Vector carrying capacity per plant k = 50,Vector birth rate 0.11.2, Vector death rate 
0.05,Disease-induced mortality rate (roguing) 0.02 [16]. 

Other plants transmission and roguing rate are 0.0064 and 0.00033 respectively|17]. Maize 
recruitment rate 276.8212,contact rate of maize plants with pathogen 8 0.0024,net decay 
rate of pathogens 0.85,contributions of I(t) to the growth of pathogen 0.0018, natural death 
rate of Sm 0.00033113 ,death rate of Im due to foliage disease 0.00099338,total population 
of 44,000 maize plants per ha and maturity of maize plants from day one the harvest day is 
151 days[8].It can also be noted that 90 days are also applicable for short season. 

Table 3 summaries the estimated and collected data as follows:- 


Table 3: Parameters values for MLND model. 














Parameters Values Units Reference/Source 
k 10 ha plant? estimated 
Bo 0.06 ha cells~'day~! [12] 

Bov 0.06 ha cells~'day~* [12] 

N, 2, 800 ha ha7! estimated 
Beat 0.286875 ha cells~'day' [13] 

Azo 0.07486 ha cells~'day~* estimated 
Ana 0.010 ha cells~'day~* estimated 
A 0.95 ha vectors™tday~! estimated 
[by 0.05 ha vectors~‘day~! 16 

Om 0.02 ha plants~'day~+ 16 

Lo 0.0005 ha pathogens 'day~! estimated 
Boni 0.0024 ha cells~'day~* 14 

Bori 0.024 ha cells~'day~* 14 

Om, 0.020 ha cells~'day~* 15 

Nm 44, 000 ha ha™! 8] 

o 0.0018 ha cells~'day~* 8] 








The numerical value of Ro referred from equation (10) is given by Ro = 7.54. 
Based on basic reproduction number in (10) and general formula of the sensitivity analysis 
in (13) we summaries the results of sensitivity indices’s obtained as tabulated in Table 4. 
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Table 4: Sensitivity indices determined based on param- 


eter values. 








Parameters Values Sensitivity indices 
k 10 -0.8122148385 
Bo 0.06 0.7996312039 

Bow 0.06 0.01258363435 
Ny 2, 800 -0.01293729282 
Dmm 0.286875 0.1869621680 
Avo 0.07486 0.01246795692 
Nii 0.010 0.7997468812 
A 0.95 0.01293729256 
[by 0.05 -0.02587458659 
Om 0.02 -0.9870627077 
lo 0.0005 -0.8122148376 
Briu 0.0024 0.00035365821 
Dori 0.024 0.00046933564 





From Table 4 we discover that the most positively sensitive parameter is the rate of interac- 
tion between contagious maize sub-population and environment Amo. In fact any movement 
of pathogens either from infected/environment to environment /infected cause much damage 
maize population. The parameter adequate or close to the most positively sensitive is the 
interaction between virus in the environment to susceptible maize sub-population 6o, to this 
juncture there is no doubt that environment plays crucial role in transmitting MLND. There 
is also positive moderate transmission of disease through maize to maize contact which is 
due to adequate contact rate of infected to susceptible maize sub-populations Bym, this is 
probably because there is a bit distance apart between one line and the other. The less 
sensitive positive parameters are A, Avo, Bov, Bmy and Bom. 

It can also be noted that, the most negatively sensitive parameter is the disease induced 
death rate om and the negatively adequate rates are natural mortality rate uo and the car- 
rying capacity of virus in the environment k, while the less negatively sensitive parameters 
are mortality natural death of vector u, and total number of vectors N». 

It is important to note that, the increase of any of positively sensitive parameter increases 
Ro and decreasing any of these decreases Ro as well, while decreasing any of the negatively 
sensitive parameter increases R, and increasing any of them decreases R, too. 

The epidemiological importance of this is to give an insights on selection of control parame- 
ters in which the most positively /negatively sensitive and adequate parameters can be chosen 
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for control strategies. 


3.6. Simulation of the basic model. In this subsection we will deal with basic model 
simulations. We will see the behavior of MLND dynamics when endemic exists and show 
how maize populations react towards vectors population and environment as well. Further- 
more, we will show the effects of most sensitive and adequate parameters in regards to basic 
reproduction number. 

All simulations are carried out using MATLAB R2013a and HP Windows 7 Professional/ 
Intel(R) Core (TM)i3-4005U CPU @1.70GHz 1.70GHz Computer machine. 


3.7. Dynamics of populations simulation. We simulate maize,vectors and viruses in the 
environment populations before looking into behavior and formation MLND based on 151 
days the maize duration from planting to the harvesting day. 

From Figures 2 and 3: We observe six different lines as an indication of various popula- 
tions. The green, cyan, black are three different lines which compose maize populations of 
at risk (susceptible), exposed and infested population respectively. The susceptible maize 
population decelerates exponentially to acquire endemic equilibrium level as they die due 
MLND which is the effects of maize to maize infection, virus in the environment and infected 
vector population. The exposed maize population assumes parabolic curve as it increases 
exponentially to a certain maximum point before exponential deceleration to the a certain 
endemic level. This behavior explain the incubation period that some resistant maize plants 
takes after exposed as it keeps increasing before turn slowly into infection class. 

The red line is the susceptible vectors and the blue line is infested vectors. The susceptible 
vectors decrease exponentially due to natural death and acquisition of infestation from se- 
verely infested maize, vectors and environment and finally acquire the endemic equilibrium 
level. The infested vectors forms parabolic curve as they do raise and drop exponentially 
to the endemicity level. The raising to the maximum point implies decrease of susceptible 
vectors add members to the infestation vectors and some might come from environment and 
infected maize. On the other hand the exponential dropping indicates no more susceptible 
vectors which means recruitment rate is much less as compared to both dying and infection 
rate. 

The virus in the environment is indicated by magenta line which as we can see extend ex- 
ponentially to a certain maximum point and proceed a constant endemic level with unusual 
number of days nearly 20 years. This behavior explain the fact that viruses in the environ- 
ment never ends even after maize extinction. The constant endemic level of viruses in the 
environment indicates that no more supply from both vector population as well as maize 
population following starvation of hunger after maize plants extinction. 
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x10 MATHEMATICAL MODEL SIMULATION VIRUSES IN THE ENVIRONMENT VS TIME(DAYS) 
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FIGURE 2. with sub-figure 2a and 2b respectively 
This Figure is showing fully simulated model and partly in the environment with the 
MLND dynamics of sub-populations/ha 


3.8. The effects of parameters on basic reproduction number. We refer Figures 4, 
5 and 6 to illustrates the effect of the most sensitive and moderate positive indexes as well 
as the most sensitive and moderate negative indexes parameters on the basic reproduction 
number. 

In Sub-figures 4a and 5b we note Amo is the most sensitive positive parameter, that the 
increase of this transmission rate result to the quick increase of the basic reproduction 
number.It imply that the interaction of infested sub-population(maize/pathogens) might 
have significant effects on maize population. Similarly, the maize death rate due to disease 
Om is the most sensitive negative parameter, that any reduction from it (stop dying of maize) 
make basic reproduction number experience significant exponential retardation. In other 
words, susceptible maize population will always remain constant if there is no infestation. 
From Sub-figures 4b and 5a we observe significant steep slope increase of basic reproduction 
number as the increase of each of these parameters (3, and Bmm as proceeds additional 
positive increments. The reason behind these rapid positive increments might be the main 
causative of MLND arises from the contact rate from infected maize to other vegetative 
plants and interaction of infested to susceptible maize rates. 

Sub-figures 6a and 6b illustrates the principal effects of natural death rate of vectors in the 
environment and the carrying capacity. Both of these moderate rate reduces exponentially 
to zero as basic reproduction number do and to the contrary is also true. This behavior 
gives an indication that killing of vectors both in the environment and its population will 
have significant effects on MLND reduction if not eradication of the whole epidemic since 
any mortality rate does not favor pathogens but maize population. 


3.9. Discussion and conclusion. The SI-SEI maize disease model which is a combination 
of both the host and vector population models has been formulated to study and analyses 
the dynamics of MLND. It is comprised of six compartments including susceptible, infected 
and exposed of maize population while of the vector is susceptible and infected with the last 


class viruses in the environment. 
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FIGURE 3. with sub-figure 3a,3b respectively 
This figure depicts the MLND dynamics partly in maize and partly in vector model 
sub-populations/ha 
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FIGURE 4. Sub-figure 4a,4b respectively 
This figure illustrates the effect of Amo and Bmm on the basic reproduction number 
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FIGURE 5. Sub-figure 5a,5b respectively 
And this show the effects of 6, and uo on the basic reproduction number R, 
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FIGURE 6. Sub-figure 6a,6b 
The last Figure which illustrates the effect of carrying capacity and maize disease induced 
death rate on the Ro 
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The disease free equilibrium (DFE) was determined which showed that in the absence of 
disease only constant susceptible maize population exists with recruitment and naturally 
dying of vectors. Other analytical results shows that the system has unique endemic equilib- 
rium point (FF) as it can be referenced from equation (12). The basic reproduction number 
Ro determined was 7.54 which show that disease is endemic, cementing the point that the 
system has unique endemic equilibrium point.From the literature review data were collected 
and summarized in Table 3 which are important for model simulations, sensitivity analysis. 
Table 4 summaries the work of sensitivity analysis with indexes that are most, adequate and 
least positively and negatively sensitive crucially for transmission of MLND. 

Finally, numerical simulation is done to illustrate the dynamics behavior of MLND in maize 
population and how some most sensitive parameters affects R,. It also suggest that viruses 
in the environment, infected maize and vector populations plays significance role in bringing 
in infection force which contribute largely to the dynamics of MLND in the maize popula- 
tion. This situation can be seen from Figures 2 and 3 where susceptible maize and vectors 
decreases exponentially while others assumes exponential curves to a certain point in time 
by increasing before they decline and find convergence point of a certain unique endemic 
equilibrium points. Another evidence is taken from Figures 4, 5 and 6 where by R, in- 
creases/decreases either with steep slopes or exponential curves as we vary by the most or 
adequate positively /negatively sensitive parameters. Lastly, the closer observation on basic 
reproduction number, sensitivity analysis and simulation indicates that Amo, Bo and Bmm 
parameters are main agents of transmission and persistence of epidemic and also very im- 
portant in finding how and where to implement the control policies for the elimination of 
MLND. The numerical solutions predict the danger of losing almost all crops if control mea- 
sures are not put into consideration since without intervention the population will vary for 
sometimes and ultimately they all keep approaching to endemic points and after a number 
of days all susceptible maize will be infected and die of disease. 

In general, we formulated and analyzed a mathematical model for dynamics of MLND for a 
single season maize populations. This model is quite simple to start with in studying and 
analyzing the dynamics of maize lethal necrosis disease in regards to maize populations as 
there are a lot of complications if we consider in depth the epidemiological process of MLND. 
Therefore study will add value to body of knowledge, policy makers and act as a bridge for 
researchers who would like to study the dynamics of MLND mathematically. 
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